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Abstract. We discuss recent numerical results obtained for gluon and ghost propagators in lattice 
Coulomb gauge and the status of the so-called Gribov-Zwanziger confinement scenario in this 
gauge. Particular emphasis will be given to the eigenvalue spectrum of the Faddeev-Popov matrix. 

Keywords: Yang-Mills theory; confinement; Green's functions; Gribov copies; Coulomb gauge 
PACS: 11.15.Hal2.38.Aw 



INTRODUCTION 

In recent years several groups have studied confinement of quarks and gluons using 
lattice simulations in Coulomb gauge [[lll-flltl. This gauge has several advantages, even 



though it breaks the Lorentz symmetry explicitly. Indeed, it is a physical gauge 111211 . 



i.e. there are no unphysical degrees of freedom. Also, a nice confinement scenario [13]- 



H19H is available for this gauge, based on Gribov 's classical work [11311 . Finally, Coulomb 
gauge is well-suited for the Hamiltonian approach and the study of hadron physics by 
variational methods 120(1 . 

Here we review lattice numerical studies in Coulomb gauge. We divide these studies 
in two periods. In the first period — called here the classical era — the studies focused 
on the infrared (IR) behavior of propagators (gluon and ghost) [111] — [0] and on the long- 
distance behavior of the color-Coulomb potential |0]-Jg]. On the other hand, in the 
second period — the modern era — the eigenvalue spectrum of the Faddeev-Popov 
(FP) operator became the main subject of investigation 



COULOMB GAUGE AT THE CLASSICAL LEVEL 



In a classical Yang-Mills theory [fl2|l . Gauss's law is written as (DjEi) a (3c,/) = p"(x,/) , 
where Z), is the gauge-covariant derivative, Ef(x,t) is the color-electric field and p^ u (x,t) 
is the quark color-charge density. (Here the sum over repeated indices is always un- 
derstood.) In Coulomb gauge, i.e. (d;A;) a (3c,/) = 0, the color-electric field can be 
decomposed into its transverse and longitudinal parts using Ef(x,t) = (E tr )f(x,t) — 
dj(j) a (x,t) , where (j) a (x,t) is the so-called color-Coulomb potential and (diE\ r ) a (x,t) = 
0. Then, Gauss's law becomes (j) a ) (x,t) = p a (x,t), where j$ = —Didi is the 
3 -dimensional FP operator and p a (x,t) = pf(x,t) - f ahc A b (x,t) (Ej r ) c (x,t) is the to- 
tal color-charge density. It follows that the color-Coulomb potential ty a (x,t) can be ex- 
pressed by means of the instantaneous and non-local operator (^~ l ) ab (x,y;t), namely 

a (3c,/) = (jZ~ l p) a (3c,/) = / d 3 y [jZ~ l ) a (x,y;t)p b (y,t) . At the same time, the 



classical Hamiltonian ffl = / <Px (Ef + Bf)/2 can be written as J4? = J#c ou i + 
Jd 3 x [(£f) 2 +^]/ 2 - Here 

^Coui = ljd 3 x (d^) 2 = \Jd 3 xJ d 3 yp a (x)y ah (xj)p b (y) , (1) 

y ab (x,y) = A)^ _1 ] a * (x,y) is the color-Coulomb-potential energy func- 

tional and we indicate with A the usual Laplacian. 

Clearly, the static color-Coulomb potential (j) a (x,t) is closely related to the 3d FP 

(equal-time) ghost propagator G(x — y,t) 8 a,b = ( Cj^~ 1 ) a (x,y;t) ) . In particular, if 

we consider the Fourier transform G(k, t) we obtain 

<TC*,0 « j d\ J d 3 y G(k,t) exp [ik } ( Xj -yj)] p a (y,t) . (2) 

Thus, if the ghost propagator has a k~ 4 = \k\~ 4 singularity at small momenta we get, in 
the limit of large separation x = |jc|, a linearly rising potential, i.e. ^"(x, t) ~ x. 



THE GRIBOV-ZWANZIGER CONFINEMENT SCENARIO 

A non-perturbative investigation of QCD in the IR limit is necessary in order to get 
an understanding of color confinement. Of course, in developing non-perturbative tech- 
niques, one has to deal with the redundant gauge degrees of freedom of the theory. The 
gauge-fixing technique developed by Faddeev and Popov assumed that one could find a 
gauge-fixing condition that uniquely determines a gauge field on each gauge orbit. How- 



ever, in Ref. 1130 Gribov showed that the Coulomb and the Landau gauge conditions do 



not fix the gauge fields uniquely, namely there existgauge-related field configurations 



that satisfy the gauge condition {Gribov copies) ||13LI19|1 . 

In order to get rid of the problem of spurious gauge copies, Gribov proposed the 
use of additional gauge conditions. In particular, for Coulomb gauge, he proposed the 
restriction of the physical configuration space (on each time-slice to the region Cl t = 
{A : (diAi) a (x,t) = 0, ^ ab (x,y;t) > 0}. Thus, inside the region Q t , the FP operator 
has no negative eigenvalues. This region is delimited by the so-called first Gribov 
horizon dQ. t , where the smallest non-trivial eigenvalue of the FP operator ^ ab (x,y;t) 
is zero. On the lattice, given a thermalized lattice configuration {U(x) } , a configuration 
belonging to the region Q. t can be obtained by finding a gauge transformation {g(x)} 
that brings the functional 4ior,f/ \s] = - L;=l T*x,t Tr u i& f ) gH x + a e h ] to 
a local minimum. Recall that, in the SU (N c ) case, both the link variables (x) and 
the gauge transformation matrices g(x) are elements of the SU (N c ) group (in the 
fundamental Af c x N c representation). 



In this review we do not discuss results related to the (possible) spontaneous symmetry breaking of the 
residual gauge freedom g(t) ||6|, 0, [IU and to the so-called A gauge i2lll . which interpolates between the 
Landau gauge (A = 1) and a Coulomb-gauge like condition (A — > 0). 



The additional gauge condition added by Gribov is not significant for the high- 
frequency vacuum fluctuations, i.e. for the perturbative regime, but it suppresses the 
low-frequency fluctuations, modifying the (non-perturbative) IR regime lfl3l [bill . In 
particular, one can show that, when the functional integration is restricted to the region 
£l t , then (on each time slice) the ghost propagator G(k,t) is IR enhanced. On the other 
hand, the transverse gluon propagator D tr (k,t) may go to zero in the IR limit, implying 
a maximal violation of reflection positivity. The latter result may be viewed as an 
indication of gluon confinement Ii22ll . Analytic results for the IR behavior of propagators 
and vertices using Dyson-Schwinger equations have been presented in Ref. Q23D . 

At the same time, the 44-component of the gluon propagator can be written |[l5|l 
as D 44 (x-%t) = V Cou i(x-y)8(t)+P(x-y,t), where V Coul (x-y) 8 ah = (V ab {%y)) 
is anti-screening and should yield a linearly rising potential, while P(x — y,t) is the 
vacuum-polarization term, i.e. it is responsible for screening and for the breaking of the 
string between color sources. One can show that these three quantities Je. a . D 44 (x — y, t), 



Vcoul( x ~ y) an d P(x — y,t)] are renormalization-group invariant [|15L ll6ll. One can also 
define the running coupling constant 

, UN c -2N f , 

S 2 Coui(k) = UN f k 2 Vcoui(k) . (3) 



c 



Clearly, if the color-Coulomb potential Vcoul (x) is linearly rising at large separation x, 
then in the IR limit we find Vc ul(k) ~ l/k 4 and gc ou [(k) ~ l/k 2 . Also, it has been 
shown llnll that the Coulomb energy of static sources is an upper bound for the static 
inter-quark potential V(x), i.e. if at large x one has Vc ul(x) — Ocouix and V(x) = ox 
then we find (N~ — l)<Jc ou i/(2N c ) > o. Analytic results for the long distance behavior 
of Vcoui(x) have been presented in Ref. II 1811 . 

Summarizing Q24D , in the Gribov-Zwanziger confinement scenario (in Coulomb 
gauge), the long-range force, responsible for color confinement, is carried by an in- 
stantaneous static color-Coulomb field. In particular, the linearly rising potential is 
related to the IR divergence of the ghost-propagator (at equal time). At the same time, 
the propagator of three-dimensionally transverse (would-be physical) gluons is ER 
suppressed and the gluons are absent from the spectrum. 



THE CLASSICAL ERA: RESULTS 

The analytic predictions described above for the gluon propagators D tr (k) and £>44(fc) 
have been verified for the SU (2) group in Refs. In particular, from Fig. 1 of Ref. 

[Q]] it is evident that, in the IR limit, the transverse propagator is suppressed, while D/^(k) 

blows up. Moreover, in the infinite- volume limit, it has been found 010] that D tr (k) is 
well described a Gribov-like propagator characterized by a pair of purely imaginary 
poles m 2 = ±iy. Numerically, at /3 = 2.2 and in the infinite-volume limit, one finds 
y = 0.33 ± 0. 14 GeV 2 . As for the ghost propagator G{k, t) , it has been studied up to now 
only in Ref. 0]. There, it is shown that G(k, t) has indeed an IR divergence stronger than 
l/k 2 . At the same time, the running coupling gc ou i(k), defined in Eq. © above, seems 



to be consistent 11,3] with an IR behavior of the type 1/k 2 . These analyses have also 
obtained o Co ui ~ o. 

In Ref. [5], the color-Coulomb potential Vc ou i(x) nas been evaluated [for the St/ (2) 
group£| as a function of the separation x, using correlators of two time-like Wilson 
lines of length 1 (in lattice units). It was found that Vc 0M /(x) increases linearly with x, 
in agreement with the 1/k 2 behavior for gc ou i(k) obtained in Refs. J3, 13]. However, in 
this case the estimate for the Coulomb string tension was Ocoui ~ 2 — 3 o. Moreover, if 
one removes the so-called center vortices [5], then the color-Coulomb potential Vc ou l( x ) 
goes to a constant at large x and Oc ou i = 0. This suggests a strong relation between these 
center vortices and the (Coulomb) confinement mechanism. Note that similar effects 
have been observed in the gluon and in the ghost propagators in Landau gauge [25] after 
removing the center vortices. 

It is also interesting that, when the temperature is turned on [@,[8j], the color-Coulomb 
potential Vqou\ (x) is not screened and it is still a linearly rising function of x. Moreover, 

the Coulomb string tension Ocoui shows a magnetic-like behavior Jsl], i.e. 0^ ul ~ 
g 2 (T)T. This implies that the Coulomb string tension cannot be used as an order 
parameter for confinement. This conclusion can be understood by observing that the 
temperature is defined by compactifying the time direction and that the Coulomb gauge 
is defined on the subspace orthogonal to the time direction. Thus, there is no reason for 
the system in Coulomb gauge to be sensitive to the deconfining transition. 



THE MODERN ERA 

In Ref. [@] the authors evaluate the gauge-field excitation energy £ (above the ground 
state energy) of a single (point-like) static color charge in Coulomb gauge. Considering 
that long-range effects should be related to the non-local-interaction term MqouU one 
finds 

£ oc y aa (x,x) = [J(- x (-ts)J(- x \ m '(x,x) . (4) 

A necessary condition for confinement is that £ should diverge in the infinite-volume 
limit, due to IR effects. (Ultraviolet divergences are regulated by the lattice cut-off.) For 
the (Coulomb) FP matrix ^ ah = —8 ab /S. — f acb A c ^ L d il one can consider (inside the Gri- 
bov region Q. t ) the eigenvalues A > and the corresponding eigenfunctions <J>^ . Then, 

Eq. © can be written as £ « (£ X F X /X 2 } with F x = V' 1 ^^)* (-A).^(4>y . 
(Here, V s is the 3d spatial volume of the lattice.) For a sufficiently large volume, 
the sums can be approximated by integrals and £ oc ( f^"' ax dX p(X)F(X)/X 2 ) , with 
/ dX p(X) = 1. In the infinite-volume limit, the volume of the Gribov region gets con- 
centrated near the Gribov horizon [0,0,122], i.e. X m i n — * 0. In the same limit, the gauge- 
field excitation energy £ blows up if 

lim"^W > o. (5) 

A^O X 



2 Similar results were obtained for the SI/ (3) group in Refs. [7, 8]. 



Thus, a necessary condition for confinement is the enhancement of p(X)F(X) at small 
momenta. 

In Appendix A of Ref. [@], an interesting analysis based on a random-matrix model 
shows that, for small eigenvalues, one should have p ( X ) = cX a if the eigenvalues X scale 

as V s 1 ^ l+a " > if the volume is increased. Numerically they find a = 0.25(5), implying 
Knin ~ l/L 2A . At the same time, they obtain F(X m i n ) ~ 1/L ~ Xf n ^ and the confinement 
criterion © is clearly fulfilled (a, llOfl. A similar result is obtained when considering 
the so-called "vortex-only" configurations [@]. On the other hand, after removing the 
center vortices, one recovers a Laplacian-like eigenvalue spectrum for the FP operator 
with lim^^Q p(X)F(X) / 'X = 0. Thus, in agreement with the findings reported in 
the previous Section, the enhancement of p(X)F(X) at small eigenvalues X and the 
confinement mechanism in Coulomb gauge seem to be strictly related to the properties of 
the center- vortex configurations. One can also show [@] that center- vortex configurations 
are (infinitely many) distinguished points on the Gribov horizon. The relation between 
these configurations and the Gribov-Zwanziger scenario in Coulomb gauge is then clear 
if, in the infinite- volume limit, the center-vortex configurations are sufficiently dense on 
the Gribov horizon. 



AN OPEN QUESTION FOR THE FUTURE 

In order to understand fully the Gribov-Zwanziger confinement scenario one should 
consider a generic gauge condition & [A] = 0, imposed by minimizing a functional E[UJ. 



Then, from the second variation of E[U], we can always define the FP matrix jtft^ . 
Clearly, when we are at a (local) minimum of E[U], the (non-trivial) eigenvalues of 
are positive and we can define a Gribov region Q. and the first Gribov horizon d£l. 
Moreover, since the configuration space has a very large dimensionality, entrop y sh ould 
favor (in the limit of large volumes) configurations near the Gribov horizon BSLliJ, 12611 . 
i.e. X m i n should go to zero in the same limit. This is indeed the case in 3d |23] and 
4d Landau gauge ll2~8ll . in 4d Coulomb gauge [@] and in 4d Maximally Abelian gauge 
(MAG) |29(]. Since the FP matrix develops a null eigenvalue at the Gribov horizon d£l, 
we should also expect the corresponding ghost propagator G(k) to blow up at small 
momenta in the infinite-volume limit. This result should in turn introduce a long-range 
effect in the theory, being probably related to the color-confinement mechanism. Indeed, 
we know from several numerical studies that the ghost propagator G{k) is IR enhanced 



in 3d [27] and 4d Landau gauge Q25|, 130, Qlfl and in 4d Coulomb gauge BQ. On the 



other hand, recent numerical results in MAG [29] suggest an IR finite G{k). Thus, the 
line of thinking reported above cannot be completely correct. Of course, one does not 
expect the ghost propagator to be particularly important for confinement in MAG, since 
in this case the accepted scenario is that confinement is related to Abelian dominance 



and (therefore) to the IR behavior of the diagonal gluon propagator [29[ 32]. In any case, 
we should try to answer the following question: what makes the ghost propagator IR 
enhanced in Coulomb and in Landau gauge but IR finite in MAG? 



A possible solution comes from the observation that in Landau 11271 12811 and in 



Coulomb gauge H] A m ,„ ~ l/L 2+a with a > 0, i.e. it goes to zero faster than in the case 



of the Laplacian. On the contrary, in MAG one has X m i n ~ 1 /L 2 ~ a (with a > 0), i.e. 
it goes to zero more slowly than for the Laplacian. This (unproven) hypothesis seems 
to be supported by the following observation. Using the same notation introduced in the 
previous section and in the limit of a large volume, we can write the ghost propagator 

G(k) as 



If we consider a FP matrix of the type ^ ab = 
Coulomb and MAG) then we have 

rbnax rlX 1 „ 



m = iJ-rLinm 2 ■ (6) 

v c a 

—S ab A — K ab (this is the case in Landau, 
$ lW = ^lI^< $ lr C7) 



In a numerical simulation we look at G(k m i n ) when the volume increases (and Xmm 
decreases). Thus, the IR behavior of G(k m i n ) depends on the quantity 

<S>°1 (k min ) = 9 1 y^ ikminX K^\ v (8) 

K min ^mm x ,y 

and there is a clear competition between the smallest eigenvalue of the Laplacian kf nin ~ 
L 2 and the smallest eigenvalue of the FP operator A„„„. 

It is interesting to notice that using Eq. © we can easily explain why finite-size 
effects are smalt] when the ghost propagator G(k) is evaluated numerically. Indeed, it 
is sufficient to have p(X)fx(k)/X ~ with /3 > — 1 in the limit of small eigenvalues 
X. In Ref. Ii28ll it has been obtained (for 4d Landau gauge) that the quantity R{X) = 

dX p(X)fx(k)X~ /G(k) behaves as X , with v > 0, for small X considering the 
two smallest nonzero momenta k. This implies j8=V — 1 > — 1. 



CONCLUSIONS 

We believe that the study of the spectral properties of the FP operator in different 
gauges can help us understand the general features of the Gribov-Zwanziger confinement 
scenario. In particular, it would be important to clarify for which gauge conditions the 
confinement mechanism can be related to an enhancement of the ghost propagator in the 
IR limit. 



3 Note that in Ref. 1I3II1 there are, actually, strong finite-size effects for the ghost propagator G(k). 
However, in that case the effects are probably due to the use of strong asymmetric lattices, with different 
ratios of the spatial over the temporal extension of the lattice. 
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